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Abstract 

We compute the high-temperature limit of the free energy for four-dimensional N = 4 
supersymmetric SU(N C ) Yang-Mills theory. At weak coupling we do so for a general 
ultrastatic background spacetime, and in the presence of slowly-varying background gauge 
fields. Using Maldacena's conjectured duality, we calculate the strong-coupling large-iV c 
expression for the special case that the three-space has constant curvature. We compare 
the two results paying particular attention to curvature corrections to the leading order 
expressions. 



1. Introduction 



There is substantial evidence to support Maldacena's conjectured duality relating 
superstring theory in anti-de Sitter (AdS) spacetime and a conformal field theory (CFT) 
[1] [2] [3] — for a comprehensive review, see ref. [4]. In particular, type IIB superstring theory 
on AdS 5 x S 5 is dual to four-dimensional N = 4 super- Yang-Mills (SYM) theory with gauge 
group SU(N C ) [1]. An interesting aspect of this duality is the study of the behavior of 
the SYM theory at finite temperature [1] [3] [5] [7] [8] [9] [10] [11] [12]. In the superstring theory 
for sufficiently high temperatures, the thermal state is described by an asymptotically 
AdS black hole [1][3]. One finds that, at a qualitative level, one recovers many of the 
expected physical properties of the thermal SYM system from the black hole geometry. A 
quantitative comparison is more difficult because supergravity provides a description of the 
SYM theory at strong coupling, and so no direct calculations can be made reliably in the 
SYM theory. In some cases [7] [14], e.g., the free energy density, perturbative calculations 
at weak coupling still reproduce the strong coupling results up to factors of order one. So 
it would appear that the corresponding coefficients are smooth functions of the effective 
coupling which interpolate between the strong and weak coupling results. In the case of 
the free energy density, subleading corrections have been calculated at both strong [10] 
and weak [11] coupling, and the results are suggestive that the interpolation may even be 
achieved by simply a monotonic function. 

Whether or not there is a smooth interpolation between strong and weak coupling is a 
question deserving close scrutiny. By carefully examining the expansions of the free energy 
at small and large coupling, and by making some simple assumptions about the behavior of 
the expansion coefficients, Li [12] has argued that it is impossible to smoothly interpolate 
between these two coupling regimes of the SYM theory. Hence Li concluded that there 
must be a phase transition at some critical coupling when the SYM theory is at any finite 
temperature. The primary motivation for the present paper was to investigate the free 
energy calculations in more detail to look for evidence of such a phase transition. As well 
as introducing a finite temperature, we consider the SYM theory on a curved background 
space. Finite temperature is included by working with a four-dimensional, ultrastatic 
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Euclidean spacetime for which the Euclidean 'time' direction is periodically identified with 
period (3 = 1/T. For the spatial geometry, we pay particular attention to three-geometries 
of constant curvature k/£ 2 (with k = 0, ±1) since these are the cases for which we may 
also compute the strong-coupling free energy using the AdS/CFT correspondence. When 
comparing the results obtained for weak and strong couplings, we follow how they depend 
on the radii, f3 and i, in the limit (3/£ <C 1 (which corresponds to the high-temperature 
limit). 

We organize the paper as follows: Section 2 states the preliminaries, describing the 
A = 4 super- Yang-Mills theory. Section 3 then addresses the weak-field calculation, which 
may be performed quite generally, using well-known heat-kernel techniques. We com- 
pute the weak coupling effective action for SYM at finite temperature, and in the pres- 
ence of slowly-varying background gauge and gravitational fields. The strong-coupling 
string/supergravity calculation is then given in section 4. A brief discussion of our results 
is given in section 5. 

2. N = 4 Super- Yang-Mills Theory 

The field content of A = 4 SYM theory [15] is {A^, \f,cp%}, where a = 1, . . . , d G runs 
over the adjoint representation of the gauge group, G, % = 1,...,4 counts the spin- half 
fields, and r = 1, . . . , 6 labels the scalars. Where necessary we will choose G = SU(N C ), 
for which d G = N% — 1, although this choice does not play an important role in the weak 
coupling calculations. 

The action of the theory may be conveniently formulated in terms of A = 1 superfields, 
of which we require one gauge multiplet, W", and three matter multiplets, 0^, all in the 
adjoint representation. The action is then given by minimal kinetic terms for all fields, 
gauge interactions, plus those interactions derived from the A = 1 superpotential: 

W=^e^C abc( p a m<p b n<p C p - (1) 

Here c a bc represent the completely antisymmetric structure constants for the gauge group. 
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This superpotential is manifestly invariant under a global SU(3) flavour symmetry 
(acting on the indices m, n and p), as well as a U(l) /^-symmetry for which the charge 
assignments of the fields are R(W2) = 1 and R(4>r) = §• These are subgroups of a 
larger SU(4) flavour symmetry which this model enjoys, defined as the automorphism 
of the supersymmetry algebra which rotates the four supersymmetry generators amongst 
themselves. Only the SU(3) x U R (1) subgroup is manifest when the theory is expressed in 
terms of N = 1 superfields. 

The theory is also scale invariant, even at the quantum level [16]. This may be argued 
within perturbation theory using the nonrenormalization theorems of N = 2 supersym- 
metry, or by constructing the most general N = 1 supersymmetric action with this field 
content which admits the above-mentioned SU(4) symmetry. 

The theory's scalar potential is 

am ^ m 

whose minima are described by the supersymmetry-preserving flat directions for which 
Cabc <Pn4>p = 0. Semiclassically, nonzero fields along these flat directions cost no energy 
and spontaneously break the gauge group, the global 577(4) symmetry and scale invariance. 
A number of Goldstone and massless gauge multiplets therefore figure prominently among 
the low-energy states, at least at weak coupling. 

In the following, we will examine the free energy of this theory at finite temperature 
and in curved background spacetimes. A nontrivial result arises because both the tem- 
perature and background curvatures typically break supersymmetry, thereby permitting a 
nonvanishing free energy. In all ways but one, the action we use for curved spacetimes is 
the same as the one just described, albeit with the substitution everywhere of covariant 
(with respect to difleomorphisms) derivatives. The only nonminimal change required is 
that the scalars couple conformally to the Ricci scalar [3]. This ensures the conformal 
invariance of the theory for general background metrics. 
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3. Weak-Coupling Calculation 



We now outline the weak-coupling calculation of the effective action and free energy 
density. Our approach is to compute the contribution of the short-distance, ultraviolet 
part of the theory to the free energy using the well-established Gilkey-DeWitt heat-kernel 
methods [19] [17] [20] [6] (see also [18]), which we will briefly review. After the heat kernel 
discussion we examine what inferences may be drawn from these methods about the long- 
distance contributions to the free energy. 

3.1) Heat-Kernel Methods 

Heat-Kernel techniques are based upon the following representation of the functional 
determinants which appear at one loop. Assuming Euclidean signature in an n-dimensional 
spacetime [19] [20] [6], 



where the upper (lower) sign is for bosons (fermions), dV is the covariant volume element, 
and K(x,y; s) satisfies the equation: 



with initial condition K(x, y; s) — > 5 n (x, y) as s — > 0. 

Part of the utility of this representation lies in the observation that the contribution 
of the short-distance part of the system to E is controlled by the small- s behaviour of 
K(x,y;s). Furthermore, this small- s behaviour has been computed once and for all, for 
general choices for the operator — □ + X. 

Concretely, if K(x, y; s) is expanded for small s: 




(3) 




(4) 



oc 




(5) 



fc=0 
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with K (x,y; s) an appropriately-chosen function [6], then E takes the following form: 

1 °° f 

where represent the following integrals: 



(6) 



Ck = f°- s k - n ' 2 exp(-m 2 s) 
Jo s 

(*"*)■ 



(7) 



= m n-2fep 



and the first few afc(x,x) are given explicitly by [17] [6]: 1 
a (x,x) = I, 

ai (x,x) = -^(R + 6X > j, 

a 2 (x, x) = -L (2 R^xpR^r - 2R^R^ + 5R 2 - 12 □ i?) 7 



(8) 



+ ^ X+ ^ 2 -^ DX + T2 y - y ^ 



In this last expression is defined by Y"^ = V M , V„ , where is the gauge and 
coordinate covariant derivative appearing in □ = g^ v V^V^. In the presence of background 
gravitational and gauge fields, this commutator becomes: 

= —iF^ v t a — — Rjju; ^ J a f3i (9) 

where t a and J a p represent the generators of gauge and Lorentz transformations on the 
field of interest. 

3.2) Applications to Spins Zero, Half and One 

We next record the above expressions for the three spins of interest in the present 
problem. For the moment we leave the dimension of spacetime arbitrary, although — with 



1 In this section, we adopt the curvature conventions of ref. [13] . 



dimensional regularization in mind — we imagine evaluating n = 4 — 2e at the end of the 
calculation. Our results are given in the presence of background metric and gauge fields. 
(Background scalars are also implicitly included through the dependence on the particle 

mass, m.) 

• Spin Zero: For spinless particles we consider the operator — □ + m\ + £ R acting on 
real scalar fields, where the choices £ = and £ = — g respectively correspond (in four 
dimensions) to a minimally-coupled and conformally-coupled scalar. Denoting the gauge 
generators as represented on scalars by t a: and assuming that we have N scalars sharing 
the same value of £ we find: 

V = ~i K» *«> ( 10 ) 

and so the a,k(x, x) become: 

tr a (x, x) = N , tr a ± (x, x) = -N ( £ + - J R, (11) 



tr o a 2 (x, x) = N 



180 



R 



\iv\p 



180 M 2 V 6 



R 2 



OR 



1 



(12) 



- -n;C(n )F"Fr. 



12 " x " u/ 

Here, and in what follows, 1Z S denotes the gauge representation of the particles of spin s, 
and the Dynkin index, C(1Z), is defined by tr(t a tb) = C(1Z) 5 a b- We normalize the gauge 
generators so that C(F) = \ for the fundamental representation of SU(N C ). 

• Spin One Half: Without loss of generality we represent spin-half particles using Majorana 
spinors. For such fields the operator of interest is^ + mi, for which the determinant is not 
well-defined. In the absence of anomalies in the Lorentz or gauge group, we instead define 
this determinant as the square root of the determinant of (^—f + mi j (jP + mi j , implying 



TYLogjV + mi) 



l - TVLog(-y 2 + m|) 
TrLog 



-□ + ml - \ R + % - Y V F« V (T alL - T* alR ) 



(13) 



Here 7 Mi , = \ [7^, 7^] and T a is the gauge generator as represented on left-handed fermions. 

The previous formalism may now be applied to this expression, provided the extra 
factor of '2' seen in eq. (13) is kept in mind. The traced commutator in this case is: 

tr {Y^Yn = ~\n, R^x P R^ Xp ~ 4C fa) F^FP, (14) 

and so the a,k(x,x) become: 

tr 1 ao(x, x) = 2iVi, tr 1 ai(x, x) = - Ni R, (15) 



tr 1 a,2(x, x) = Ni 



Again Ni denotes the total number of Majorana spinors involved. 



(16) 



• Spin One: Gauge potentials, A a , are the fields representing spin-one particles. Their 
contribution to the one-loop effective action takes a form for which the above-described for- 
malism applies provided that the gauge is chosen appropriately. For background-covariant 
Feynman gauge the appropriate differential operator for A® has the form — 5^ 5% □ + 
8% {m\)l + X££ , with X^ and the traced commutator, Y"^, given in this case by: 

Xg = -SSR£ + 2iF<* v (r c ) a b 

(17) 

tr (Y^Y"") = -N t R^ Xp R^ Xp + (A) F^Ff . 

Here (r a ) c = ic b ac denotes the gauge generator in the adjoint representation, and C(A) is 
the corresponding Dynkin index. (For SU(N C ) we have C(A) = N c .) 

In addition to the vector potentials, each spin-one particle also requires a Fadeev- 
Popov-DeWitt ghost, which contributes like a minimally-coupled complex scalar in the 
adjoint representation, however it is an anticommuting field. Combining the vector and 
ghost contributions gives, for Ni spin-one particles: 

2 

tr 1 ao(x, x) = 2N±, tr 1 a±(x, x) = — N\ R, (18) 



tr i ci2(x, x) = N\ 



13 22 5 I 

— R,ux P R^ Xp + ^ R^Br --R 2 +-UR 



(19) 



3.3) Combining Results 

We now proceed to the case of interest, N = 4 SYM theory in four spacetime dimen- 
sions. This corresponds to the choice Nq = 6iV~i, Ni = 4N\ with N± = d G given by the 
dimension of the gauge group. We also choose all particles to transform in the adjoint rep- 
resentation, TZq = IZi = A, and we take £ = ~, as is appropriate for conformally-coupled 
scalars. For this choice, and combining according to Tr = tr — tr i au + tr i a^, we 
find: 

Tr a (x,x) = JV - 2M + 2JVi = Tr ai (x, x) = -6 d G [ £ + - J R = 0, (20) 



Tr a 2 (x,a;) = d G ^ R^R^ + 



3U + - 



# 2 



2 



1 



ti^vti — -K 



(21) 



Note that Tr a 2 takes the form familiar from the conformal anomaly in which the coefficient 
oiR^ Xp R^ Xp vanishes [22] 2 



3.4) UV Finiteness on Ultrastatic Spacetimes 



As may be seen from the form of the integrals, Ck, of eqs. (6) and (7), the coefficients 
Tr dk(x,x), for k = 0,1,2 determine the ultraviolet divergences of the four-dimensional 
theory (n = 4). The vanishing of the k = and k = 1 terms in eqn. (20) and the fact 
that the expression for the k = 2 term in eqn. (21) depends only on the background 

We thank Arkady Tseytlin for pointing out an error in our original version of this formula. 
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curvature (due to the cancellation among different spins of terms proportional to F^ U F^ U ) 
is a check on the calculation. It shows that this theory, on flat spacetimes, is ultraviolet 
finite at one loop, as expected. We now show that this cancellation also occurs for arbitrary 
four-dimensional ultrastatic spacetimes. 

The finite-temperature calculation of the next section is performed for ultrastatic 
spacetimes. These admit metrics for which coordinates may be locally chosen to ensure 
ds 2 = dr 2 +"f mn dx m dx n , where the spatial metric, 7 mn , is independent of r. Since any such 
metric is effectively three-dimensional, both its squared Weyl tensor and its Gauss-Bonnet 
integrand must vanish identically: 



Eq. (22) may be used to eliminate the quantities R pv \ p R puXp and R^R^ in terms of R 2 , 
with the result: 



It is now obvious that eq. (21) vanishes when eq. (23) are used. Thus we see that 
N = A SYM is UV finite at one loop on all four-dimensional ultrastatic spacetimes. This 
result also extends to background geometries which have a product structure, e.g., S 2 xS 2 . 
Again, since the individual components of the metric are lower dimensional in such a case, 
one finds the vanishing of eq. (22). However, note that for more general backgrounds, it is 
quite possible for Tr a,2(x,x) to be nonzero. For instance, for four-dimensional maximally 
symmetric spaces R^\ p = R (g^\9u P - g PP 9v\), and so AR^R^ = R 2 , giving 



e P P e Pi R pva pR\p^g 



R^x P R^ Xp - 2 R^R^ + l -R 2 = 
-R^x P R^ Xp + 4R^R^ -R 2 = 0. 



(22) 




3 



(23) 



Tva 2 (x,x) = -^R 2 ^0. 



(24) 



In such cases logarithmic UV divergences would appear in the effective action. 



3.5) Nonzero Temperature 
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The calculation of the ultraviolet part of the one-loop effective action for nonzero 
temperature is performed along lines very similar to the zero-temperature result just de- 
scribed [21] [18]. We do so by restricting to the ultrastatic metrics of the previous section, 
and periodically identifying the Euclidean time r with period (3 = 1/T. The only change 
required in the previous calculation is to find the Heat Kernel which is periodic — or, for 
fermions, antiperiodic — in time. 

Such a kernel may be obtained from the one used above at zero temperature through 
the method of images: 



K±(t-t',x,x.';s)= (±) n K(t-t' + n/3,x,x';s). 



(25) 



where, again, the upper (lower) sign is for bosons (fermions). This sum may be evaluated 
within the small-s expansion for the ultrastatic spacetimes of interest here because the 
time-dependence of Kq(t — t',x, x; s) is quite simple: 



K (t - r', x, x; s) = K (xl, x; s) exp 



(r-r') 



t\2l 



As 



(26) 



Expanding K± for small s, as in eq. (5), and using the known form for Kg, one finds [21]: 



K ± (x,x;s)=K (^x;s) £ (±)»expf-^-J 

T) — — nun \ * 



(27) 



The sums may be performed, giving standard Jacobi t?-functions [23] : 

K + (x,x;s) = # 3 (r) K (x., x;s), K-(x, x; s) = ^ 4 (r) K (x., x; s), (28) 
where r = i(3 2 /(4ns). 

Substituting eq. (28) into expression (3) for the effective action, one then finds that 
the free energy density is: 



1 

F = T W77Z W2 J2 C i ±} tr a kM, 



2 (47r) n A 



(29) 



fc=0 
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where now CY represent the following integrals: 




k-n/2 



k-n/2 



exp(— m 2 s) t?3 



exp(— m 2 s) $4 




(30) 



These integrals diverge for k = 0,1 and 2 if n = 4 due to the limit of integration 
at s — > 0. This infinity is an ultraviolet divergence, and is removed if one focuses purely 
on the temperature-dependent part of the problem. This is accomplished by computing 
A^ — cj^ — Ck, which gives the change of the free energy relative to the zero-temperature 
result: AF = F(T) - F(0). 

The potential infrared divergence as s — > oo does not arise so long as m is nonzero, 
so A^ may be evaluated numerically to any desired accuracy. For our later purposes, 
however, it is instructive to examine the high-temperature limit, m<T. This limit can 
be subtle due to the shadow cast by incipient infrared divergences which can re-emerge in 
this regime. A well-known example of this occurs in flat space, since Aq" 1 "-* acquires terms 
such as m 3 T when expanded in powers of m/T [24], [25]. This nonanalytic dependence of 
Aq + ^ on m 2 as m 2 — * thwarts its evaluation via term-by-term integration after expanding 
its integrand in powers of m 2 , since the successive integrals which are obtained in this way 
diverge more and more severely in the infrared. Such nonanalytic dependence on m 2 is 
characteristic of the singularities which arise in the presence of massless particles. 

For nontrivial background fields we see that evaluating at m = causes the A^ to 
diverge at s — > oo for sufficiently large k. This is very similar to the divergences which 
are found when a series in powers of m 2 is attempted in the absence of background fields. 
In this case the divergences are a consequence of a breakdown of the derivative expansion 
due to the nonanalytic dependence on the background fields which is generated as the 
curvatures cut off the infrared singularities which are generated by the various massless 
modes [27]. 

In four dimensions these divergences potentially arise at m = when k > 2. We 
now explicitly display the first of these by evaluating the terms k < 2 of the Gilkey-De 
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Witt expansion. One finds a result expressed in terms of Euler's T-function and Riemann's 
^-function: 



d s k-n/2 



8 

oo 



Jo s 



0s 



4tts 
4tvs 



4 \ 2 K 



2 y <-*)«„-»), 



_ ( l-n+2k 



(31) 



= (2 



1) A« 



Combining all expressions, and expanding about n = 4, gives the following result for 
the free energy density of a general theory: 



AF = 



1 



[tr B a k (x,x) + (l-2 1 - n+2k ) tr F a k (x,x)], 



2 (4tt)^/ 2 



fc=0 



7T 2 T 4 

90 
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tr B a + - tr F a 
2 



^2 

24 



32tt^ 



2 In ( — ] + 

4-n VT 0y ' 



tr B ai +- tr F ai 



1 31n2 

[tr B a 2 - tr F a 2 J + — T tr F a 2 + 



(32) 



Here tr B a k and tr F respectively denote the trace over the bosons and fermions of the 
model. The ellipses within the square brackets represent terms which neither diverge, nor 
involve the logarithm of T. By contrast, the ellipses at the end of the equation represent 
terms which involve higher than two powers of curvature and gauge field strengths. 

The infrared divergence arises in eq. (32) as the pole as n — > 4 in the k = 2 term 
of the series. Notice that the condition for the cancelling of this particular divergence is 
tr B ci2 = tr F a,2, which is precisely the condition for the cancelling of the zero temperature 
logarithmic ultraviolet divergence. As a consequence this particular divergence does not 
arise in the present example of N = 4 SYM, although the same is not expected to be true 
for the infrared divergences arising for higher k. 

Evaluating this expression using the a k s of N = 4 SYM, and using the simplifications, 
eq. (23), which follow from the restriction to an ultrastatic spacetime, finally gives: 



AF = 



n 2 d G T 4 
6 



24 



1607T 2 

13 



8tt 2 



[iv 1 - a 



(33) 



In the next section we compute the free energy density in the strong-coupling limit 
using the AdS/CFT correspondence [1]. We do so for the special case of vanishing back- 
ground gauge fields, and for the specific case where the spatial three-geometry has constant 
curvature The Ricci scalar for such a space may be written as R = —6k/£ 2 , where £ is 
the radius of curvature, and k = +1, 0, —1 for a three-sphere, flat space and a hyperbolic 
three-plane. Further for comparison purposes, we will be interested in the gauge group 
for SU(N C ) in the large-iV c limit, so d G = N 2 — 1 N 2 . For these choices, our result in 
eq. (33) yields: 

7T 2 T 4 KT 2 



AF = Nt 



6 

tt 2 N 2 T 4 
6 



A£ 2 



+ o 



3k 
2n 2 £ 2 T 2 



£^T 2 ) 



+ o 



1 



(34) 



4. Strong-Coupling Calculation 



Using the AdS/CFT correspondence, the dual supergravity description of the N = 4 
SYM theory at finite temperature is an asymptotically anti-de Sitter black hole [1] [5] . An 
appropriate class of metrics describing Euclideanized black holes is [28] 



(35) 



where k = +1, 0, —1 and dT,s(k) is the line element on a three-dimensional manifold with 
constant curvature k/£ 2 . An explicit representation of the latter may be chosen as 



gE 3 (k) 



(1 - Kp 2 )dz 2 + 



dp 2 



1 — Kp' 



+ p 2 d(P 2 



(36) 



but the precise form of the three-dimensional metric will not be needed in the following. All 
of these five-dimensional metrics (35) satisfy the equation of motionf R^ v = —(4/£ 2 )g fIiy , 



Note that there is a change of conventions for the curvatures between this section and the previous one. 
Those of this section follow those of ref. [26], since these are conventional in the gravity literature. Although 
the metric signature remains — +++, we must replace R^vXp^ — R^xp in all formulae. 
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which arises from the action 



/ = 



16ttG 

To avoid a conical singularity at 



^- j d 5 x^(R + 12/ t 2 ) - -L- J d 4 xVhK . (37) 



r+ = ^ (j/ K 2 £ 2 + 4r%e 2 - K^j , 



(38) 



where Qtt vanishes, one must select the period of r to be 



1 

Now, the Euclidean action evaluated for the classical black hole metric is interpreted as 
giving the leading contribution to the free energy. In the context of AdS black holes, such 
calculations were first carried out in ref. [29] . As typically arises in these calculations, how- 
ever, this action diverges and so to produce a finite result, we subtract off the contribution 
of a reference metric [30]. This step is analogous to subtracting off the zero-temperature 
free energy in the field theory calculation — see section 3.5. In this case, the reference 
geometry is simply anti-de Sitter space which is produced by setting ro = in eq. (35). 
Taking care to match the asymptotic geometries correctly, one finds 

A/ = ^AF=^^_±A. (40 ) 

Note that in defining AI = V S /3AF : the relevant volume is essentially the coordinate 
volume of r and dE^K). It is not the proper volume of a surface of constant radius at 
large r in eq. (35), which of course diverges as r — > oo. 

Now the duality prescription [1] identifies G 5 £ 5 = 87r 3 g 2 a /4: and £ 4 = 4irgN c a' 2 , 
where g and a' are the string coupling constant and the inverse string tension, respectively. 
Combining these formula then yields the free energy density at strong coupling for SU(N C ) 
SYM in the large- N c limit: 



8 V i 2 
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where the function J-{x) is given by: 



1-3 



16 



1 




(42) 



where the final expansion applies for x « 1. This expansion corresponds precisely to the 
high temperature limit in eq. (41), which then becomes 



The latter is now in a form which is readily compared to the weak coupling result (34) . 

For the k = case, of course, the leading term is the full answer up to corrections 
in the effective coupling. For the k = +1 and -1 cases, similar calculations yielding the 
curvature corrections appear in refs. [14] and [31], respectively. 

Notice that for k = or -1, AI in eq. (40) is always negative, and so the black hole 
solution (35) always provides the dominant saddle-point in the supergravity path integral 
for any temperature. On the other hand, for k = +1, AI becomes positive for r+ < £, 
which corresponds to ft > 2tt£/3. Hence in this low temperature range, AdS space itself 
is the dominant saddle-point. This change at ft = 2n£/3 is believed to be associated with 
a deconfinement phase transition in the SYM theory [5]. 4 For comparision to the weak- 
coupling calculation, we are interested in the high temperature regime where the black 
hole is the dominant AdS saddle-point for all three choices, k = 0, ±1. 



In comparing the free energy densities calculated at one loop (34) and strong coupling 
(41), we see that these expressions, while not identical, yield two expansions which look 

4 Note that below this phase transition at P=tt£/V2, there is a branch point in JF(/3 2 /I 2 ). This corresponds 
to the minimum temperature for which a black hole solution with k=+1 exists. 




(43) 



5. Discussion 
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very similar. Both give a free energy having a high-temperature limit which has the form 
of a Taylor expansion in powers of f3 2 /£ 2 : 



where the coefficients b n (X) are functions of the effective 't Hooft coupling, A = g\ M N c = 



The significance of this observation for the strong-coupling expansion is in what 
eq. (44) does not contain. In particular, the weak-coupling expansion had the poten- 
tial to involve both fractional powers of (3 2 /£ 2 (similar to the m s T term in Aq + ' ) ), as well 
as terms proportional to log/3 (such as when Tra 2 7^ 0). Such terms reflect the infrared 
divergences which are associated with the massless bosons of the perturbative spectrum. 
They do not arise at one loop in the weak-coupling calculation due to the cancellations 
in Tra2, which reflect the finiteness and conformal invariance of N = 4 SYM theory in 
ultrastatic background geometries. This same cancellation does not hold in general for 
all terms in the weak-coupling limit. For instance, even on flat space the massless gauge 
bosons generate infrared divergences to perturbative calculations of thermal quantities like 
gluon damping rates [32], once one proceeds beyond one loop. The same might also be 
expected to follow for the k > 3 terms of the derivative expansion at one loop. 

Unlike the weak-coupling result, the strong-coupling expression, eq. (42) is not simply 
given as an asymptotic form for large temperatures, but is explicitly given in terms of the 
function JF(k/3 2 /£ 2 ) for the geometries of interest. Explicit examination of J-{x) clearly 
shows it to be analytic in its argument near x = 0, indicating the absence of singularities 
in the high-temperature regime. But the existence of exact conformal invariance of the 
underlying model ensures the absence of a gap in the particle spectrum. In four dimensions 
this follows directly if the ground state respects the conformal invariance, or from the 
existence of a Goldstone 'dilaton' mode if the conformal symmetry is spontaneously broken. 
Infrared divergences would be expected to be weakened or absent if the only massless 
states were Goldstone modes [33], since these decouple in the long- wavelength limit [34]. 
The absence of infrared singularities in the large temperature limit of the strong-coupling 




(44) 



/(2«' 2 ). 
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calculation therefore indicates either the existence of only Goldstone massless modes, or 
the persistence into the large- iV c limit of the one-loop cancellations of infrared divergences 
which were seen at weak coupling for non-Goldstone massless modes. 

Notice there is no singularity predicted by eq. (42) for any other temperatures, in 
spite of the branch point which apparently arises at x = n 2 /2. This branch point has no 
direct implications because: (i) x = —(3 2 /£ 2 < if k = — 1; (ii) the function J-(x) collapses 
to a constant (and hence the branch point vanishes) if k = 0; and (Hi) the branch point 
lies outside the high temperature domain of validity ((3 < 2n£/3) of the strong-coupling 
calculation if k = +1. 

The behaviour we find is consistent with the behaviour which was previously argued 
for N = 4 SYM by Witten [5]. It would be of interest to extend other techniques which 
have proven powerful in flat space, such as those of ref. [35] , to analyse the curved spaces 
considered here in more detail. In particular, so far as they go our results are consistent 
with a QCD-like picture in which no phase transition occurs for any finite coupling, A c . 

Although our calculation cannot reveal the functional form for b n (X) away from the 
limits A — * and A — > oo, it does differ in some ways from similar calculations of the 
free energy density in the Higgs phase of the SYM theory [8]. In that case, the results 
suggest that the strong-to-weak coupling interpolation is achieved with a single overall 
multiplicative function. While one might have hoped for a similar simple form in the 
present calculation, it is clear from our expressions that this is not the case. Rather, the 
weak coupling form (34) determines 

6 (0) = 1, 6i(0) =| , 6 2 (0) = 0, (45) 

while the strong coupling expansion (41) fixes 

^(A-oo)-^, ft^A-oo)-^, 6 2 (A^oo)^. (46) 

Following the arguments presented in ref. [14], the free energy densities calculated 
using Euclidean techniques above also have an interpretation as Casimir energy densities. 

18 



The latter arises if the Wick rotation back to a Minkowski signature manifold is made on 
an appropriate 'time' coordinate in the constant curvature three-manifold. Then eqs. (34) 
and (41) yield the Casimir energy density for SYM on dSs x S , R 3 x S 1 and AdS% x 5" 1 for 
k = +1, and —1, respectively. The negative Casimir energy is generated by antiperiodic 
boundary conditions on the 5" 1 factor for the fermions. 

Finally, in closing, we remark that the calculations in section 3 show that in a general 
curved space, N = 4 super- Yang-Mills need not be finite. We observed though that for 
background geometries with a product structure into two (or more) lower dimensional 
components, UV finiteness is maintained due to a remarkable cancellation of terms in 
Tra2, which is given in eq. (21). Indeed, in the finite temperature calculations, finiteness 
arose because of the ultrastatic form of the background metric, which gives a product 
structure of the form S 1 x M3. Of course, it is not difficult to find backgrounds where no 
such cancellation occurs. For example, consider the simple case of S 4 , which is maximally 
symmetric with R = 12 /£ 2 where £ is the radius of curvature. In this case eq. (24) implies 
the nonzero result 



and hence a logarithmic UV divergence appears in the SYM effective action. 

It is not difficult to choose coordinates on five-dimensional AdS space such that the 
'boundary' manifold (i.e., the asymptotic regulator surfaces) is S 4 , i.e., 



where dfl 4 is the standard metric on a unit four-sphere. Now the logarithmic UV di- 
vergence appearing at one-loop in the SYM effective action would at first sight seem to 
present a problem for the proposed AdS/CFT correspondence, but in fact, it provides a 
remarkable consistency check [37]. The AdS/CFT correspondence provides a new set of 
intrinsic surface terms [36], [37] (for related work, see [38]) for the gravitational action (37). 
Remarkably, these new terms render the gravitational action finite except precisely in the 
case where Tra2 7^ for the boundary surface. In the latter case, the volume integral 




(47) 



ds 2 = j + r 2 dQ 4 

1 4- — 



(48) 
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produces a logarithmic divergence instead. For the S 4 boundary, one can precisely match 
this divergence with that appearing in the dual field theory calculation [37]. Hence, in 
keeping with the UV/IR relation [39] of the AdS/CFT correspondence, one again finds 
that an ultraviolet divergence in the field theory is matched by an infrared divergence in 
the supergravity. 
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